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The following report serves as the technical progress report for 
Contract NAS9-15476 which is entitled, "Analysis of Scanner Data for 
Crop Inventories" . This report describes the work carried out under a 
set of SIX tasks for the period 15 March 1978 through 14 June 1978. 

Work on this contract is performed in the Infrared & Optics Divi- 
sion directed by Mr. Richard R. Legault. Dr. Qunntin A. Holmes, Head 
of the Information Systems and Analysis Department is the Project 
Director for thxs contract and Mr. Richard F. Nalepka, Head of the 
Multispectral Analysis Section, is the Principal Investigator. 

This contract is part of a comprehensive and continuing program 
of research at ERIM into earth resources information systems which 
employ remote sensing of the environment from aircraft and satellites 
to gather data and which use automatic processing to extract infor- 
mation from the data. The basic objective of this multidisciplinary 
program is to develop such information systems and practical tools 
which will provide planners and decision-makers extensive accurate 
information quickly and economically. 

The six tasks reported on for this reporting period include: 

Multisegment Training 

Evaluation of Partitions for Signature Extension 
Wheat Vs. Small Grains 

Forecasting Production of ^flieat from Satellite Data 
Prepilot Study of Multicrop Spectral Separability 
Multicrop Labeling Aids 

No report is Included for the Haze and Soil Correction or the 
Evaluation of Multitemporal Classifiers tasks since, for this quarter, 
the sponsor has redirected our efforts to multicrop associated activi- 
ties. 
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TASK 1 

MULTISEGMENT TRAINING 
(R. KautL and W. Richardson) 


1.1 OBJECTIVE 

The objective of Task 1 is to develop a sampling strategy for 
selecting trainxng data, applicable to proportion estimation over a 
wide region. The main requisites of that strategy are that it produces 
a representative sample and that the training sample size is small com- 
pared to the total area to be classified. 

1.2 APPROACH 

1. Create a conceptual basis for the problem of training in a 
large scale remote sensing system, incorporating the inputs 
from UCB, LARS, and other ERIM tasks, and consistent with 
LACIE operational constraints. 

2. Within this framework, propose a detailed methodology for 
training selection. 

3. Demonstrate the selection methodology in an intermediate scale 
exercise over a partition containing from 15 to 30 sample seg- 
ments from which 5 to 10 segments are selected for training and 
for xjhich a \?heat proportion estimate is made. 

4. Incorporate both multitemporal and across partition signature 
extension capability into the final procedure. 

5. Incorporate the capability to work with incomplete sets of 
multi temporal data and to optimize selection to make estimates 
at several times during the growing season. 
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1.3 PROGRESS 

Durxng the last quarter a baseline version of Procedure Bhyas 
provided to Task 2 of this contract for testing. The approach to the 
missing acquisxtion capahxlity was further clarified and coding modi- 
fications needed were identifxed. Coding of the missing acquisition 
capahilxty is about 90% complete. The problem of defining a composite 
procedure combining desirable aspects of both Procedure B and JSC’s 
Procedure 1 was considered. Ro final design was defined but some 
suggestions were included in the SRS.T quarterly review, June 12 - 16. 

ERIM personnel took part in the Procedure 2 design review, June 16, 1978. 

The major effort during the quarter was in the development and 
exercising of diagnostic tools and procedures to measure the performance 
of components of Procedure B. The tools developed are similar to some 
already developed at JSC to measure component performance of Procedure 1. 
Note that the tests planned on the baseline version of Procedure B under 
Task 2 are tests of global performance compared to other approaches. 

The component performance tests being developed and carried out under 
this task are for the purpose of identifying and isolating the sources 
of variance and bias in Procedure B and of establishing optimal para- 
meter settings . 

1.4 TECHNICAL DISCUSSION 

The technical discussion will be limited to the development and 
exercise of the component performance tests. Some introductory discussion 
is followed by the material presented at the quarterly reviex/. 

1.4.1 COMPONENTS AND PERFORblANCE MEASURES 

Currentljr, three major components of Procedure B are being examined: 

1. Spectral/Ancillary Data Stratification (B-CLUSTERING) 

2. Training Segment Selection 

3. Training Blob Selection 
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The performance measures which are being used to evaluate these 
components are the variance, the bias, and the number of training seg- 
ments required to represent the entire set of segments. 

1.4.2 SPEGTRAL/ANCILLARY DATA STRATIFICATION 

The main performance measures to be used in evaluating the perform- 
ance of B-CLUSTERING are the reduction in variance of the proportion 

i 

estimate due to the clustering and the number of training segments 
required to achieve 90% of maximum "value" in segment selection. 

The formula for reduction of variance due to stratification is 


f = 


Fn.p. (1 - p,) 
11 1 

NpCl - p) 


where the sum is over all B-clusters 

til 

is the number of pixels in the i B- cluster and N is 

the total number, 

is the proportion of \rheat in the i^^ B-cluster and p is 
the average proportion of wheat, 

f is a number which ranges from 0 to 1, and is better if 
it is smaller. 

The reduction of variance is computed using ground truth sites 
so that the quantities and p^ are kno^m. 

The reduction in variance of the B-CLUSTERIKG is affected by the 
spectral channels used, the ancillary variables used, and in particular 
by the relative weights used. As the weight on the ancillary variables 
is increased, the stratification becomes better but the number of seg- 
ments required to represent the entire partition increases. Hence the 
best choice of ancillary and spectral variables is the one with the 
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smallest f value which achieves a given percent of representation with 
a certain fraction of the segments used for training. 

1.4.3 TRAINING SEGMENT SELECTION 

The performance measures to he used in evaluating the training 
segment selection procedure are the variance and bias of the partition- 
wide proportion estimate and the individual segment proportion estimates. 
The parameters which are varied are the number of segments chosen for 
training and the number of blobs chosen for training. In particular, 
i-t is of interest to know whether the process of selecting segments 
intrinsically Introduces a bias into the procedure. If so, it may be 
necessary to randomly select samples from all segments i Even if this 
turns out to be the case. Procedure B may have sampling advantages over 
a single segment procedure. 

1.4.4 TRAINING BLOB SELECTION 

The performance measures for training blob selection are again 
the variance and bias of partition-wide proportion estimates and of 
individual segment estimates. The parameters which can be varied are 
the minimum number of internal pixels and the shape factor (compactness) 
in blobs chosen for training. 

In order to measure variance and bias, it is necessary to create 
random replicates of the selection process. This has been done by 
making different random choices of the blobs actually used in the 
procedure. 

1.4.5 QUARTERLY REVIEW PRESENTATION 

The following presentation includes recommendations and future 
plans . 
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OUTLINE 


Introduction 

Status 

Experimental Results 

Conclusions 

Recommendations 



INTRODUCTION TO PROCEDURE B 


Procedure B is a technique of proportion estimation which tells an analyst 

WHICH scene elements TO LABEL^ AND THEN USES THOSE LABELS IN AN UNBIASED WAY 
TO PRODUCE A PROPORTION ESTIMATE FOR THE SCENE 

Fundamental Concept Similar to Procedure 1 

- Stratified sampling used for "bias correction" 

" Sampling approximately proportional to size of strata 
minimizes variance 

- Strata formed by "classification" of samples using a 
clustering algorithm operating on a feature vector 

Major Differences in Concept 

- Preprocessing/feature extraction 

- Number of strata 

- Number of segments 



PREPROCESSING/FEATURE EXTRACTION 


« Data Screening 
® Sun Angle and Haze Correction 
® Features Used 

- Spectral: Brightness/Greenness 

~ Temporal: 2 or 3 Biowindows Simultaneously 

- Spatial: Contiguous Groups of Spectrally-Temporally 

Similar Pixels (BLOBS) 

- Ancillary Data: 

Parameters such as view angle> crop calendar^ 
available soil moisture^ latitude and LONGITUDE/ 
WHICH vary from SEGMENT TO SEGMENT 


TASK 1 
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STRATIFICATION PROCEDURE 


® Cluster Together Defined Spectral-Spatial-Temporal Features 

- Multi segment 1/2 State of Kansas) 

- Ancillary Data (Values Specific to Each Segment) 

- Multiple Strata (B-clusters) 


TASK 1 



TRAINING SELECTION PROCEDURE 


Procedure has evolved to an allocation procedure for stratified 

SAMPLING 

Attempt to select a subset of segments which will provide a sufficient 

SUPPLY of samples FOR LABELING WITHIN EACH STRATA 
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PROPORTION ESTIRIATION PROCEDURE 


® Estimated proportion is the weighted average of the estimated 

PROPORTIONS IN ALL THE B-CLUSTERSj WHERE THE WEIGHTS ARE THE 
NUMBER OF PIXELS IN EACH B-CLUSTER 

* The estimated proportion in each B-CLUSTER is the weighted average 
OF the estimated proportions of each blob used for labeling^ where 
the weights are the numbers of pixels in each blob 

® The estimated proportions in each BLOB come from analysts (ground 
truth) in production (research) version 


TASK 1 
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PROCEDURE B PROPORTION ESTIMATION FORMULA 


A 

P 



I 


WHERE 

Nj IS THE TOTAL NUMBER OF PIXELS IN THE lIH bCLUSTER 


A 





k^t(i) 


^IK “"iK 


z: 


B 


K£T(i) 


IK 


AND WHERE 

t(i) is THE SET OF BLOB CHOSEN FOR TRA-ININ6 THE I—- BCLUSTER 

IS THE NUMBER OF PIXELS IN THE I BLOB USED FOR TRAINING 

Pjj^ IS THE ESTIMATED PROPORTION FOR THE IK™ BLOB 


TASK 1 
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PROCEDURE B DATA FLOW 


Feature 

Extraction 


Clustering 


Training 

Selection 


Proportion 

Estimation 



TASK 1 
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STATUS OF PROCEDURE B DEVELOPMENT 


° HAVE ESTABLISHED A BASELINE PROCEDURE FOR T&E 

° HAVE INITIATED EXPERIMENTS TO MEASURE PERFORMANCE 
OF THE COMPONENTS OF PROCEDURE B 


TASK I 



EXPERIMENTS 


Objective 

- DETERMINE OPTIMAL PARAMETER SETTINGS 
FOR COMPONENTS OF PROCEDURE B 

Approach 

" MEASURE component PERFORMANCE AS FUNCTION 
OF PARAMETER VALUES 

Measures of component performance 

- VARIANCE reduction FACTOR Cr.V.) 

- SAMPLING BIAS 


SAMPLING VARIANCE 



COMPONENTS BEING EXAMINED 


spectral/spatial stratification (blob) 
spectral/ancillary data stratification CB~Cluster) 
training segment selection 
training blob selection 


proportion estimation algorithm 
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SELECTION/PROPORTION ESTIHATION 

“ Parameters to vary or effects to evaluate 
" Segment selection 

: NUMBER OF SEGMENTS CHOSEN 

: VARIATION DUE TO CHOICE 

- BLOB SELECTION 

: Number of BLOBS chosen 

: Variation due to choice 

- Proportion estimation 

: Algorithm detail 

» Performance measures are bias and variance of 
proportion estimate for the entire group 
OF segments 


TASK 1 
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WHEAT PROPORTION (PERCENT) 



9 6^1 

number of segments used for training 

Procedure B: Proportion estimates for 9 segments with approximately 300 blobs 

SELECTED FOR TRAINING AND WITH NO ANCILLARY VARIABLE. 
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WHEAT PROPORTION (peRCENT) 



NUMBER OF SEGMENTS USED FOR TRAINING 

Procedure B: Proportion estimates for 9 segments with approximately 600 blobs 

SELECTED FOR TRAINING AND WITH NO ANCILLARY VARIABLE, 
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NUMBER OF SEGMENTS USED FOR TRAINING 


Procedure Bj Proportion estimates for 9 segments with approximately 1500 blobs 

SELECTED FOR TRAINING AND WITH NO ANCILLARY VARIABLE. 



WHEAT PROPORTION (PERCENT) 



Procedure B; Proportion estimates for 9 segments with approximately 3000 blobs 

SELECTED FOR TRAINING AND WITH NO ANCILLARY VARIABLE. 



V'/HEAT. PROPORTION (.PERCENT) 



NUMBER OF SEGMENTS USED FOR TRAINING 

Procedure B: Proportion estimates for 9 segments with approximately 300 blobs 

SELECTED FOR TRAINING AND WITH NOVEMBER SOIL MOISTURE, 
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Procedure B: Proportion estimates for 9 segments with approximately 600 blobs 

SELECTED FOR TRAINING AND WITH ilOVEMBER SOIL MOISTURE. 



WHEAT PROPORTION (.PERCENT) 
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Procedure B; Proportion estimates for 9 segments with approximately 3000 blobs 

SELECTED ?"OR TRAINING WITH NOVEMBER SOIL MOISTURE, 
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NUMBER OF SEGMENTS USED FOR TRAINING 


Procedure B: Proportion estimates for 9 segments with approximately 1500 blobs 

SELECTED FOR TRAINING AND WITH NOVEMBER SOIL MOISTURE. 



SPECTRAL/ANCILLARY DATA STRATIFICATION 


) • 
ro 

CO 


» Parameters to be varied 

- SPECTRAL VARIABLES USED 

- ANCILLARY VARIABLES USED 

- RELATIVE WEIGHTING 


^ Performance measure is reduction of variance due to stratification 


F 


XjNi Pi (1-P^) 


N P(l-P) 
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Proportion of Hheat Bias^ (P - P). 
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RELATIVE WEIGHTING ON 6 SPECTRAL VARIABLES C GREENESS 
AND BRIGHTNESS FOR 3 BIOWINDOWS) 

- VARY ALL 6 VALUES 

" NORMALIZE TO 90 BCLUSTERS 

“ METHOD OF STEEPEST DESCENT 

" 2 WAY R.V. 

- FOUND BEST OPERATING POINT WHEN VARIABLES 
WERE WEIGHTED INVERSELY TO THEIR RANGE' 

- 1 WAY R.V, = .539 
~ 2 WAY R.V, = .406 


TASK 1 
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3 GREEN VARIABLES COMPARED WITH 3 GREEN AND 3 BRIGHTNESS VARIABLES 


1 36 + 3B 

36 

I WAY 

R.V. 

1 

.622 

2 WAY 

R.V. 

.406 

.430 
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ONE-WAY R.V. 




Reduction of Variance factor due to clustering 

NOVEMBER 1 SOIL MOISTURE 

ALONG WITH THE SPECTRAL VARIABLES, 

T MEASURES THE WEIGHT PUT ON THE ANCILLARY 
VARIABLE, '^ = 72 CORRESPONDS TO ZERO WEIGHT. 


70 


80 


90 100 
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NUMBER OF 
CLUSTERS 
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Reduction of Variance factor due 'to 

CLUSTERING 

6 SPECTRAL VARIABLES 

AS A FUNCTION OF THE NUMBER OF CLUSTERS 












Segj^ent INDE)^ number 


0 180 200 


80 100 120 140 160 180 200 


Variance Reduction Factor as a Function of Number of BCLUSTERS FOR SINGLE SEGMENTS’.' Variance 

REDUCTION FACTOR DUE TO BLOBBING IS SHOWN AS A LOWER BOUND. SLIDE 1 OF 2 , 
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VARIANCE REDUCTION FACTOR. F. VS MAXIMUM MISREGISTRATION 

ESTIMATE (PIXELS) 
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VARIANCE REDUCTION FACTOR, F, vs flUfIBER OF. BLOBS 
WHICH SURVIVE AFTER STRIPPING 


Correlation .,7 
SiSNIFICANCE .04 





2pi 


FORMEHUY WILLOW RUN LABORATORIES THE UNIVERSITY OF MICHIOAN 


TASK 2 

EVALUATION OF PARTITIONS FOR SIGNATURE EXTENSION 
(R. C. Cicone)* 


2.1 INTRODUCTION 

The sampling and classification strategy of the Large Area Crop 

Inventory Experiment (LACIE) entails employing local signature training 

to determine wheat proportion estimates within 5x6-mile sample segments, 

Wheat proportions are then aggregated within designated strata. Multi- 

/ 

segment signature extension is philosophically founded on the premise 
that representative training information may be determined using non- 
local procedures at an additional savings in cost and reduction in the 
variance of the estimate. 

Task 2 is concerned with addressing the key issues found in 
Table 2.1 that pertain to non-local training techniques. 

2.2 OBJECTIVE 

The objective of Task 2 is to test and evaluate techniques and 
procedures which embody the signature extension approach to large area 
crop inventories using Landsat data. 

2.3 APPROACH 

The approach adopted to address the objective of Task 2 is out- 
lined in Table 2.2 

2.4 SUMMARY OF PROGRESS DURING SECOND QUARTER 

Table 2.3 reviews progress of Task 2 during the first quarter of 
the contract. Progress during the second quarter is outlined in 
Table 2.4. 

A 

T. Wessling and J. Stinson contributed to the work reported. 
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2.5 DETAILS OF PROGRESS 

Efforts this quarter have concentrated on a parametric evaluation 
of multisegment signature extension in North Dakota using a technique 
called preclassification. Figure 2.1 diagrams the effort carried out. 
Results will be described for cases employing Biowindows 1 and 2 in 
this spring wheat growing state, uncorrected and sun angle/haze corrected 
data, without partitioning of segments. 

The preclassification technique employed is briefly described in 
Table 2.5. A more detailed description of this technique is available 
in Reference 1. The North Dakota data set used in this analysis is 
illustrated in Figure 2.2. 

Of primary interest to this task is the role of partitioning in 
a multisegment signature extension environment. Figures 2.3(a), (b) , 
and (c) illustrate a partitioning of North Dakota according to ten-year 
averages of precipitation and degree day parameters [2] . The analysis 
initiated in '.this quarter is concerned with techniques used to make 
use of static partitions. Several approaches are outlined in Table 2.6. 

Figures 2.4 and 2.5 illustrate analysis conducted to date in an 
unstratified multisegment signature extension environment. 

The ordinate axis of Figure 2.4 represents classification accu- 
racies; the ordinate axis of Figure 2.5 represents bias in 209 dot 
grain proportion estimates. The abscissa specifies the number of seg- 
ments used for training. Computations were carried out using" different 
combinations of training segments drawn randomly T/ithout replacement. 

The graphs contain the average result with one standard deviation bar 
drawn about the average. 

Figures 2.4(a) and 2.4(b) illustrate field center classification 
accuracy for grain and non^grain classes. Figure 2.4(a) was derived 
from uncorrected Landsat data. Figure 2.4(b) illustrates an improved 
grain classification accuracy using haze and sun angle corrected data 
[3]. Note the overall poor average classification accuracy achieved 
for. grains mthout stratification of the segments. 
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Figures 2.5(a) and 2.5(b) illustrate the bias in 209 dot grain 
proportion estimate. Thfe variance indicated is that due to the random 
draw of training segments, not that due to segment-to-segment variance 
in accuracy. The uncorrected case, (Figure 2.5(a)) illustrates a tend- 
ency to underestimate grain by about 1% in the first case to over 7% 
in the last. The bias in the XSTAR corrected data is more stable 
across all cases beginning at about 1.7% to just over 3%. In any case 
the classifier displayed high variance segment-to-segment accuracies 
(not illustrated) . 

Corrected data seems to be establishing a pattern of improved 
results. This is supported also by the analysis reported last quarter 
using Kansas field mean data. It is yet to be seen what improvements 
can be made in a stratified environment. It is clear, however, that 
in an unstratified environment, multisegment- signature extension is 
at a loss in coping with overlap of spectral distribution of different 
real classes that is most lihely attributable to a number of ancillary 
conditions . 

2.6 CONCLUSIONS AND RECOMMENDATIONS 

Tables 2.7 and 2.8 contain conclusions and recommendations drawn 
based on the analyses conducted through the first two quarters of this 
task. 


2.7 PLANS 

Table 2.9 lists the activities to be conducted through the 
remainder of this contract year. 
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TABLE 2.1 


KEY ISSUES 


“ Multisegment Area Estimation Procedures 
“Signature Extension 
“ Use of Static Partition Boundaries 
» Data Normalization and Preprocessing 


TASK 2 



TABLE 2.2 


APPROACH 

Test Design and Implementation Considerations 

- Experiment design 

- Program development 

- Data preparation 

Parametric Evaluation Using Preclassification 

“ Training gain 

- Haz-e correction 

- Static partition boundaries 

- Training segment selection 

Procedure Evaluation 

- Local Procedure B 
■ - Procedure B 

- Multisegment Procedure 1 



TABLE 2.3 


REVIEW OF PREVIOUS QUARTER'S PROGRESS 

o Evaluated RIultisegment Signature Extension 

- Using 1975-76 Kansas field means data base 

- Preclassification technique 

®XSTAR Haze and Sun Angle Correction Resulted in Substantial Improvement 
in Field Center Classification Accuracy 

©The Use of UCB Static Degree Days and Precipitation Strata Significantly 
Improved Results^ at a Cost in Training Gain 

©The Use of a Green Development Indicator as a Segment (Hatching Criterion 
in Biowindow 1 Resulted in Improved Classification Accuracy 



Task 2 
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TABLE 2.4 


SU.WRY OF PROGRESS 


® Program Development 

- PRECLAS. TASCAP, STRIP. CO^PRS. GNDTRll. ILEC 

® Prepared North Dakota Data Set 

“ 15 Lacie Phase HI Blind Sites 
4 Biowindows 

- JSC Ground Truth 

o E valuation of I^ultisegment Signature Extension Using Preclassification 

- Proportion Estimation + Field Center Accuracy 

- 2 and 3 Biowindows (15 and 11 Segments) 

- Corrected and Raw Data 

- UCB Partitions 


Task 2 
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TABLE 2.5 


PRECLASSIFICATION 

© A Sum of Likelihoods Classifier 

“ Cost-Effective Means of Testing 

- Can Examine Significance of a priori Heighting Factors 

- Facilitates Operational Use of Partitions with 
Hultisegment Procedures 

©Each Segment Classified by Each Training Segment Separately Using Sum 
of Likelihoods, Quadratic Classifier 

©Decision Made by Comparing Appropriately Weighted Sums of Wheat and 
Other Likelihoods from Different Training Segments 


Task 2 
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TABLE 2.6 


USE OF PARTITIONS 


o 


UCB Degree 


Days^ Precipitation 


©Strategies Considered 

- Treat boundaries as static 

- Weight training from same strata more heavily 
“ Limit training to same or adjacent strata 


- Use partition variables as distance measures for weighting training 

- Insure each partition is represented by a training segment 

oUsE Function of Distance in PIiles Between Segments for'Weighting Signatures 



Task 2 


553 
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TABLE 2.7 


CONCLUSIONS 


° Unstratafied signature extension is affected by undesirable 

"SEGMENT-TO-SEGMEMT OVERLAP BETWEEN GRAIN AND OTHER SIGNATURES^' 
WITH CLASSIFICATION DOMINATED BY MORE PREVALENT OTHER CLASSES. 


* Preprocessing techniques appear to be un-doing the partial overlap 
OF otherwise separable GRAIN AND OTHER PIXELS WHICH IS DUE TO 
SCENE-TO-SCENE haze AND SUN ANGLE VARIABILITY, 
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TABLE 2.8 


RECOmENDATIONS 


® In two multi segment tests conducted even the best results in an 

UNSTRATIFIED ENVIRONMENT WERE POOR, An UNDERSTANDING OF HOW TO 
USE STRATAFICATION TO IMPROVE THESE RESULTS MUST FIRST BE ATTAINED 
BEFORE USING MULTISEGMENT SIGNATURE EXTENSION OPERATIONALLY, 

® Haze and sun angle preprocessing should be incorporated into any 

PROCEDURE EMPLOYING SIGNATURE EXTENSION TECHNIQUES, 
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TABLE 2.9 


PLANS 

° Examine Role of Static Stratification 
» Document Results 

» Test Local Procedure B and Compare to 
Procedure 1 Approach 

“Examine Uses of AflOEBA in Procedure B 

“Test Multisegment Procedure B 
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FIGURE 2.1 


EVALUATION OF MUIJISEG^ENT SIGNATURE EXTENSION 

PROCEDURE 
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FIGURE 2.2 


NORTH DAKOTA DATA SET 
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figure 2.3(b) 
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figure 2.3(c) 
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FIGimE 2.5(a) 


MULTISEGFIEHT SIGNATURE EXTENSION 
15 N.D. Phase III Blind Sites 
Biowindows 1 & 2 

209 Dot Proportion Estimation Accuracy 
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FIGURE 2.3(b) 


XSTAR DATA 
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TASK 4 

SPECTRAL SEPARABILITY OF SPRING TTOEAT FROM OTHER SMALL GRAINS 
(¥,A. MAIILA AND J.M. GLEASON*) 


4.1 INTRODUCTION 

The problem of distinguishing bet^jeen spring wheat and other small 
grains is a problem of interest to LACIE and similar agricultural crop 
inventories. A prior study at ERIM** gave indications of separability 
of spring wheat and barley under certain conditions in Landsat multispectral 
data from several Phase 2 LACIE Blind Sites in North Dakota. Last 
quarter, we reported confirmation and explanation of many observed 
characteristics in the Landsat signatures of small grains, through 
analysis of 1975-76 LACIE Field Measurement data and USDA crop reports 
and crop production statistica. Also several implications for discrimination 
procedures were dratra from analysis of the results. 

The objective of this Task is to develop a spectral classification 
method for discriminating spring wheat from other small grains, using 
Landsat data. 

4.2 APPROACH 

Table 4.1 summarizes the approach taken toward meeting the 
Task objective. 

4.3 SUPtMARY OF PROGRESS DURING THE SECOND QUARTER 

During the reporting period, two data sets from the 1976-77 
growing season v^ere prepared and analyzed. First were LACIE field 
measurements of the reflectance of spring small grains, made in 
Williston, N. Dakota, using an Exotech Model 100 (Landsat-band) 


* ¥. Malila served as Task leader; E.P. Crist also contributed to the 
reported work. 

** Malila, William A. and James M. Gleason, "Investigations of Spectral 
Separability of Small Grains, Early Season WTieat Detection, and Multi- 
crop Inventory Planning, "ERIPf 122700-34-F, Environmental Research 
Institute of Michigan, Ann Arbor, Michigan, November 1977. 
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radiometer. 'Second t?ere LACIE Phase 3 Blind Site data from North and 
South Dakota, Montana, and Minnesota. Substantial data preparation 
activities were conducted to produce a comprehensive Phase 3 data set - 
included were screening, corrections for haze and sun angle, sampling 
on a 5x5 grid (897 points), and merging \7ith crop labels extracted from 
JSC ground truth tapes. 

Trends observed in LACIE Phase 2 data and corresponding field > 
measurements were also found in these more recent data sets. Unitempora’’- 
correct classification percentages for spring wheat vs. barley were in 
the 80 's given complete training on individual segments and 76% for a 
multisepment example. Spring x^heat and barley x^ere most separable in 
the grox'7th stage x^hen they turn color from green to yelloxj-green or yellox^ 
(the soft dough stage) . Barley tends to turn color sooner than xJieat and 
some of it has a distinctive brightness after turning is complete. 

Tasseled - cap brightness again X'7as the key feature for discrimination; 
together with greenness it captured most data variation and spectral 
differences. A parameter indicating crop maturity v7ould appear promising 
for exploring differences in crop calendar , xjhere coupled x^ith crop calendar 
estimates. 

The best linear discriminant tended to overestimate the ratio of 
x-7pring x-7heat to barley, even given complete training on the small grain 
observations. 

Finally, from this our first analysis of extensive Landsat data 
on oats, X7e found that the spectral separation of sprang x^heat from oats 
to be substantially more difficult than from barley. 
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4.4 DETAILS OF SECOND QUARTER ACTIVITIES 

4.4.1 ANALYSIS OF FIELD MEASUREMENT DATA 

The data set analyzed is described in Table 4.2. It represents a 
portion of the LACIE Field Measurements made .in North Dakota during the 
1976-77 growing season. The instrument has filters which cover the 
spectral bands of Landsat. Linear discriminant analysis was performed 
on these measurements as a function of observation date. The results 
presented in Table 4.3 represent an upper bound on separability that 
might be achieved. Note that separability was best after heading. 

The next part of our analysis was to examine LACIE Phase 3 Landsat 
data directly. 

4.4.2 ANALYSIS OF LACIE PHASE 3 LANDSAT DATA 

An extensive data set ^?as assembled by processing 'data provided by 
NASA/JSC. Landsat data from 28 segments in four states (See Figure 4.1) 
were prepared in a manner described in Table 4.4. Preprocessing included 
corrections for atmospheric haze and sun angle and a systematic sampling 
of the pixels. The 897 points selected include the 209 dots of Procedure 
1 as a s'ubset. 

Twelve segments with substantial numbers of barley pixels were 
analyzed for (linear discriminant) separability of labeled spring wheat 
and barley using Landsat spectral data. The results presented in 
Table 4.5 for individual segments by acquisition date show clearly 
that in 1977 the July 12-13 time period (Julian Day 193-194) was the 
best date for separability. Most correct classification values are in 
the 80% range. Each number in the table is the conditional average 
correct classification, i.e., the average of the spring v?heat and barley 
correct classification values. It was noted that the spring wheat value 
generally was larger than the barley value. 
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We wished to rank order 'the tasseled~cap features in their con- 
tribution to separating spring wheat and bailey. Seven segments 
(excluding Montana 1929) had acquisitions on the best date. First, 
their features were rank ordered on an individual segment basis, with 
results as shown in Table 4.6.. Next, haze-corrected data values were 
pooled, a composite multisegraent linear decision rule was established, 
and the tasseled cap features were ranked. Brightness was most important 
in both the composite and most individual segments. The second- rankdd 
feature varied from segment to segment, but greenness was selected for 
the composite data set. 

The location of the multisepment decision line in Greenness-Brightness 
space IS shom in Figure 4.2 x^ith respect to all spring wheat points and 
in Figure 4.3 with respect to all barley points. The classification 
accuracies achieved by this line (four- dimensional version) are 
presented on the right half of Table 4.6. Overall, x^heat was 85% correct, 
barley 68%, and the average 76% correct. Values for individual segments 
also are given. Note the atypical results obtained for segment 1699 and 
the explanation in the footnote. This result highlights the need for a 
good crop calendar estimate to help establish decision Dines betx^^een 
spring wheat and barley. 

In addition, xi?e examined the spectral separability of oats and spring 
x^heat in Tx^elve segments. (Our LACIE Phase 2 data set did not have 
sufficient oats for earlier comparisons.) The results presented in Table 4.7 
indicate both poorer separability than spring xjheat vs. barley and no 
clear best date for separation. 

Finally, xxre performed multitemporal clustering of spring wheat and 
barley points for one segment (1663) to better examine temporal trajectories 
and their variability. Plots for four xvheat cDusters and txTO barD.ey clusters 
are presented in Figures 4.4 and 4.5, respectively. 
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4 . 5 •RECO^MENDATION 

It is recommended that crop growth stage information be included xn 
the information acquired as part of the periodic observations of fields 
in blind sites. 


4 . 6 PLANS 

Flans for the third quarter are presented in Table 4.8. 
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TABLE 4.1 


APPROACH 

o Gain Understanding of Potential Discrimination Features, Such As Differing 
Spectral Characteristics Vs. Crop Calendar and Differing Crop Calendars 

" By using available Landsat data 

- Supplemented by other data 

o© LAC IE field measurement data 

®®USDA agricultural statistics 

® Establish, Test and Evaluate One or More Classification Procedures 

® If Results Sufficiently Promising, Establish Sensitivity of Classi- 
fication Performance to Acquisition Sequence and Accuracy of Training 

® Begin Consideration of Advanced Procedures, e.g,, Year-to-Year or use of 
Thematic Mapper Bands 


TASK A 
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TABLE 4-2 


SHALL GRAINS REFLECTANCE DATA 

1977 LACIE Field Measurements in Williams County^ North Dakota^ 

With Exotech Model 100' Radiometer 


o 32 Plots 

A Crops: Hard Red Spring Wheat 

Durum Spring Wheat 
Barley 
Oats 

2 Varieties Each Crop 

2 Soil Moisture Levels: Fallow in 1976 

Wheat in 1976 

’2 Replicates Each Combination 
® 6 Measurement Dates 


Spring Wheat 



Julian 

Growth Stage 

m 

-< 

to Growth Stage 

6/18 

169 

2 

2 

Tillering 

6/23 

17A 


3 

Stem Extension 

7/03 

18A 

3. A 

A 

Heading and Flowering 

7/lA 

195 

A 

5 

Ripening 

7/20 

201 

5. A 



8/05 

217 

5 
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TABLE 4.3 


LINEAR DISCRIMINATION RESULTS ON 1977 RADIOMETER REFLECTANCE DATA 
SPRING WHEAT (HARD RED AND DURUM) VS. OTHER SMALL GRAINS (BARLEY AND OATS) 


Decision 

Class 


JuN 18 (Day 169) 

True Class 


sw 


1 

03 


SW 


OSG 


1 

55 1 

1 

23 

A5 , 

77 


OSG 

No, Observations 29 30 


AVG. 

CORRECT 


66 


Classification Percentage 

JuN 23 (Day 17A) 

True Class 


sw 

OSG 


sw 


OSG 


75 

25 

25 

75 


16 16 


AVG, 

CORRECT 


75 


JuL 3 (Day ISA) 
True Class 


sw 

OSG 


sw 


OSG 


81 

31 

19 

69 


16 16 


AVG. 

CORRECT 


75 


Decision 

Class 


Jul 1A (Day 195) 
. True Class 


sw 

OSG 


sw 


No. Observations 58 


OSG 


83 

13 

17 

87 


60 


AVG. 

CORRECT 


85 


Jul 20 (Day 201) 
True Class 


SW 

OSG 


sw 


9A 

6 

6 

9A 


OSG 


16 16 


AVG, 

CORRECT 


9A 


Aug 5 (Day 217: 
True Class 


sw 


OSG 



3/A OF Oats and 1/2 Barley Harvested by Aug. 5, 
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TABLE 4.4 


SPRING SMALL GRAINS DATA SET PREPARED FROM LACIE PHASE 3 BLIND SITES (1976-77) 

® 28 Segments From Four States: 16 North Dakota 

6 South Dakota 
A Montana 
2 Minnesota 

® Preprocessed Landsat Data 

- Screened visually and digitally to flag clouds^ shadows^ 

AND bad data 

“ Acquisitions were selected 

- Corrected for haze and sun angle (XSTAR Algorithm) 

- Transformed to Tasselled-Cap Features 

® Sampled 897 Points (5x5 Grid^ Includes 209 Dots as Subset) 

o Merged With Ground Truth Data From JSC Tapes (Labeled as either 
a single specific cover code or a mixture pixel) 

0 Associated Ancillary Data 

- LACIE crop calendar estimates 

- Latitude and longitude 
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TABLE 4.5 


SPRING WHEAT VS. BARLEY DISCRIMINATION 

Landsat Data from LACIE Phase 3 Blind Sites (from 897 points) 
Complete training using ground truth.- 
Linear decision rule. 


Average Percent Correct on Date; 




No. Pixels 

May* 

May 

JUN 

JUN 

JUL 

JuL 

- Aug 

State 

Segment 


B 

1-3 

18-21 

5-8 

23-25 12-13 

29-31 

17-18 

Minn 

1515 

287 142 



63 

56 

85 

** 





1523 

139 

75 

59 

66 

64 

56 



68 


67 

ND 

1616 

288 185 

59 

• 67 

63 






62 


1619 

315 

67 

55 

55 

67 

72 



60 


65 


1622 

252 103 

60 

66 

61 

57 





56 


1637 

194 

28 

65 

55 

54 


83 





1640 

224 103 

60 

69 


64 

78 


78 


69 


1663 

239 

93 

60 

64 

73 

69 

77 


65 


55 


1899 

308 224 

52 

61 

65 

62 1 







1927 

127 

47 

57 

67 

64 

71 1 

^51 




58 

SD 

1699 

54 

23 


72 

66 

89 1 





64 

Mont 

1929* 

61 

21 

72 

70 


70 

69 


88 




Landsat acquisitions over Montana were approximately one week after other acquisitions', 


Box denotes largest value for segment, 
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TABLE 4.6 


RANKING OF TASSELLEL CAP FEATURES FOR DISCRIFlINATING BETWEEN 

SPRING WHEAT AND BARLEY 

July 12,13 Acquisitions (1973) 

(Days 193,194) 

XSTAR-Corrected Landsat Data 


Segment 

Order of 
Choice 

Crop 

Calendar 

MULTISEGMENT RULE 

% Correct Classification 
Sp Wht Barley AVg. 

^INN 1515 

B,6,Y,N 


5.4 

90 

77 

84 

SD 1699 

B,N,G,Y 


6.0 

o«- 

100 

50 

ND 1637 

N,B,G,Y 

• 

5.3 

92 

54 

73 

.1640 

G,B,Y,N 


5.4 

88 

61 

75 

1663 

B,6,Y,N 


5.3 

83 

74 

78 

1899 

B,G,N,Y 


5.3 

93 

58 

75 

1927 

B,G,Y,N 


5.4 

81 

81 

81 

All Seven I 

b.gj.yI 



85 

68 

76 


Note advanced crop calendar for Segment 1699 which caused spring wheat 
TO MOVE TO BARLEY SIDE OF DaY 193 DECISION LINE; On DaY 176 WITH CROP 
CALENDAR 5,1, THE RESULT WAS SW = 96%, B = 48%, AND AVG, = 72%, 
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TABLE 4 


SPRING WHEAT VS. 




Landsat 

Data from 

LAC IE 

Phase 3 



Complete 

training 

USING 

GROUND ■ 



Linear decision rule. 




No. Pixels 

May 

May 

State 

Segment 


B 

1-5 

18-23 

Minn 

1515 

287 

36 




1523 

139 

58 

55 

56 

ND 

1602 

215 

22 

55 

55 


1637 

194 

22 

67 

61 


1652 

116 

21 

55 

64 


1661 

188 

38 

52 



1663 

239 

27 

54 

66 


1903 

63 

27 

62 



1927 

127 

37 

62 

64 

SD 

1675 

85 

21 

59 

64 


1686 

27 

84 


IZSJ 


1699 

5^ 

66 


72 
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TABLE 4.8 

PLANS 

“ Conduct More Extensive Studies of the Relationship Between Periodic 
Ground Observations and Landsat Data for Wheat in Phase 3 Blind S'ites 

® Continue Analysis of Field Measurement Data 

^ Develop and Test a Spring Wheat Vs. Barley Discrimination Procedure 
Based on a Spectral Crop Maturity Index 
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FIGURE 4.1 


LOCATIONS OF LACIE PHASE 3 BLIND SITES ANALYZED 
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FIGURE 4.2 

SCATTER OF SPRING WHEAT DATA 
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SCATTER OF BARLEY DATA 
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BRIGHTNESS 



TEMPORAL TRACKS OF S 


WHEAT CLUSTER TRAJECTORIES 5 DATES 
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FIGURE 4.4 


TE.WRAL TRACKS OF SPRING WHEAT CLUSTERS (SEGMENT 1663) (Cont.) 
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TASK 6 

FORECASTING PRODUCTION OF TOEAT FROM SATELLITE DATA 

(R.F. Nalepka) 


6.1 OBJECTIVE 

The objective of this task is to investigate the utility of 
Landsat, meteorological, and ancillary data for forecasting winter 
wheat yield. 

6.2 APPROACH 

The approach to fulfilling the objective of this task is to examine 
two separate topics related to the general xssue of wheat production r 
forecastine. The principal topic to be investigated is the construction 
and implementation of improved and combination yield models. Another 
topic to be investigated is the generality and large area applicability 
of various yield forecasting approaches. In order to satisfactorily 
investigate the two topics we are also developing a larger data base 
with which to test various approaches. 

6.3 TECHNICAL DISCUSSION 

Our current' technical activities are described in Figure 6.1. 
Details of the generation of new data are sho\ra in Figures 6.2 and 6.3. 
The preliminary data analysis that has been done is shown in Figure 6.4. 

Figure 6.5 shows a comparison of Landsat, CCEA, and early KCLRS 
forecasts of yield for counties in the southwest Crop Reporting District 
(CRD). KCLRS estimates of 1976 yield made in March of 1977 are used as 
the standard of comparison. 

Figure 6.6 shows a stepwise regression analysis of the relative 
utility of various factors for estimating yield in 13 counties of 
southwest Kansas. 
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The results xnvolving Landsat estimates of green measure and 
yield reported here should be considered as preliminary in nature. 

The optimal procedure to get an accurate estimate of Landsat green 
measure and yield is not yet clear, and there are indications of some 
error (although compensating) in the methods we used to generate this 
data. A more controlled method would he to use data from knox'm xjheat 
fields, on blind sites. However, it is not yet clear that such a 
procedure \<rill produce a representative sample of xjheat conditions. 

It has been proposed that one factor limits crop yield in any 
particular situation, and that this concept could provide the basis for 
a yield model that is superior in some respects to multiple linear 
regression. A law— of-minimum (LOM) model can be constructed by deter- 
mining relationships betx^een separate factors and yield under conditions 
identified as limiting to each factor. The model can be implemented by 
calculating predicted yields from all of the above relationships, and 
choosing the minimum value for each observation as the correct prediction 
of yield. 

We have constructed LOM models on data from 13 counties in southxi?est 
Kansas using two different procedures, and compared these models with 
a model generated using multiple linear regression. The results are 
shown in Figure 6.7. 

These results suggest that an LOM model may be better or worse than 
multiple linear regression, depending on exactly hox^ the model xjas 
constructed. Our preliminary conclusion is that construction of an 
optimum LOM model probably can't be consistently achieved using a 
single unsupervised algorithm. Rather, considerable human interaction, 
judgement, and perhaps trial and error, may be required. 
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Figure 6.1 

Current*4ctivities 

CTi 

I 

• Generate Data Set of Landsat^ Meteorological and Cultural 
Variables for Southwest Kansas, similar to Data Set Previously 
Established for Central CRD. 

«* Preliminary analysis of New Data (including Cate*- Liebig) 
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Figure 6.2 

GENERATION OF LANDSAT DATA 


Make estimates of Landsat Green Measure on 

COUNTY'-BY'-COUNTY BASIS 

Make Landsat estimates of yield on county-by-county 
BASIS USING Finney ITS data for training 



Figure 6.3 


GENERATION OF ANCILLARY DATA SET 

Organize Monthly Meteorological Data 

Compute CCEA yield estimates 

Calculate KCLRS yield estimates 

Determine Cultural Data (% summer fallow, % irrigated) 



Figure 6.4 


PRELIMINARY ANALYSIS OF DATA 

Compare Landsat and CCEA estimates of Yield with KCLRS Yield 

Generate law-of-Minimum models for counties of Southwest Kansas 

USING Landsat Green Measure^ Aug - Feb precipitation^ May Precipitation^ 

percent summer fallow^ percent irrigated 

Generate Linear Regression Model using above data 

Compare LOM and Regression Models 
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Figure 6.5 


COMPARISON OF SEPARATE LANDSAT AND CCEA YIELD ESTIMATES FOR 
COUNTIES IN THE SOUTHWEST CRD 
(compared to March 1977 KCLRS estimates). 


Esxjmaxe Mem 

CCEA ‘ 30.0 

(1 July 1976) • 

Lands AT 28.5 

(18 April 1976) 

KCLRS 17.0 

(April 1976) 


%AN Diff from KCLRS 
+3.0 


+1,5 


- 10.0 
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Figure 6.6 


REUTIVE UTILITY OF VARIOUS FACTORS FOR ESTIMATING YIELD 
(13 COUNTIES OF SW Kansas assessed by Stepwise Regression) 


STEP 

VARIABLE 

r2 

S.E. 
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.62 
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Figure 6 . 7 

COMBINATION YIELD MODELS 
(based on 13 Counties in $W Kansas) 

Model .S-iLl 

Multiple Linear Regression 0,76 1.6 

LaW“OF~Minimun (Fit R) 0.65 1.7 

LaW“OF-Minimun (Fit M) 0.86 1,0 
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TASK 7 

PREPILOT STUDY OF MULTICROP SPECTRAL SEPARABILITY 
( W.A. Malila and E.P. Crist)* 

7.1 INTRODUCTION 

Through a re-organization of efforts, making them better focussed on 
LACIE transition objectives, this task was formed and is being reported 
here for the first time. The key issue being addressed is the extension 
of large area crop inventory technology to include important crops in 
addition to wheat. 

7.2 OBJECTIVE 

The objective is to conduct prepilot signature studies using currently 
available data, to provide insights and identify potential problem areas 
for investigation in future multxcrop pilot testing and P^.OID activities by 
JSC and SR&T contractors. 

7.3 APPROACH 

The approach being taken toward meeting the objective is out lined 
in Table 7.1. Note the emphasis on use of currently available data. 

7.4 SUMMARY OF PROGRESS DURING THE QUARTER* 

Table 7.2 summarizes our progress during the reporting period. 

7.5 DETAILS OF PROGRESS 
7.5.1 ANALYSIS OF CITARS DATA 

Table 7,3 describes the CITARS data set which was the primary focus 
of our first quarter’s activity. The locations of the sites are shown in 
Figure 7.1. As indicated in Table 7.4, both field-center pixels and Landsat 
mean vectors formed from them are being analyzed. 

As our first step in investigating corn and. soybeans signature 
characteristics, we conducted a series of unitemporal and multitemporal 


* W. Malila is task leader; contributions to the reported work were 
also make by J.M. Gleason. 
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clustering operations on four-date data from the Fayette, Illinois, segment. 
The results are presented next in the fonn of spectral signature plots 
of Landsat Band 6 vs. Band 5. The dispersion of pixel values within each 
cluster is indicated by a "one— sigma" ellipse about the mean value. 

The first series of plots is of multitemporal (four-date) clusters for 
corn arid soybeans. All clusters with more than five pixels are illustrated. 
Figures 7.2 arid 7.3 show the way corn pixels progress from a June 10th 
' dispersion along a line, indicating much visible bare soil (of varied 
reflectance), to compact distributions on July 17th and August 21st. In 
contrast, the soybeans clusters shown in Figures 7.4 and 7.5 are widely 
dispersed on all dates, with most compactness on the last date, August 21. 

A contribating factor is the fact that some soybeans were planted several 
weeks late in fields from which winter wheat had been harvested. 

The patterns for these crops are shown in more detail on a uniteraporal 
basis in Figures 7.6 through 7.9; all clusters with ten or more pixels 
are shown. Clustering for the latter two figures included pixels from other 
crops as well as corn and soybeans. 

Finally, three-date multitemporal clusters, generated using all labeled 
pixels, are presented in Figures 7.10 and 7.11. 

From these initial results, we do not attempt to draw conclusions 
regarding the best date for classifying corn and soybeans or the accuracy 
levels to be expected. However, we do note that substantial spectral 
differences do exist, that the distribution of corn on the last two dates 
is compact, and that multitemporal clusters were more pure than unitemporal 
clusters. 

In addition to the preceding Landsat data analysis, we have digitized 
the periodic (18-day) observation data from the Fayette segment, to enable 
subsequent comparisons of spectral data and agronomic characteristics of 
the fields. 
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7.5.2 ANALYSIS OF LACIE PHASE 3 BLIND SITE DATA 

To increase the data base for corn and soybeans, we have prepared 
data from eleven LACIE Phase 3 Blind Sites which contain some corn or 
soybean fields. The data preparation steps and amounts of data 
available are summarized in Table 7.5. 

7.6 PLANS 

Plans for the next quarter are presented in Table 7.6. 
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TABLE 7.1 


APPROACH 

• Analyze Available Landsat Data on Corn and Soybeans 

- 1973 CITARS Data (Illinois and Indiana) 

Re-examine in light of understanding and new techniques 
gained through LACIEi 

•• Structure of Landsat data (Tasselled Cap Transformation) 
•• Spectral/spatial and multitemporal clustering 
••XSTAR preprocessing (corrections for sun angle and haze) 
•• Quantitative use of periodic ground observations 

- 1977 LACIE Blind Site Data (Selected Sites on Fringes of a 
U.S. Corn Belt) 

• Analyze Available Data from Other Sources 

- Field Measurements of Spectral Reflectance 

- Simulation Model Data 
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TABLE 7.2 


SUMMARY OF PROGRESS 

• Have Begun Re-Examination of CITARS Data in light of recent technology 
- Landsat Signatures as Function of Time 

MULf'lTMPORAL CLUSTERING OF FaYETTE SEGMENT 
“ Quantitative Analysis of Periodic Ground Observations 


® Have Completed Preparation of Eleven Phase 3 Blind Sites Containing 
Some Corn and Soybeans for Multi crop Analysis 

'• Have Requested Earlier Field Measurements of Corn and Soybean Spectral 
Reflectance Data from LARS 



TABLE 7.3 


Cl TARS DATA SET 


Landsat Data 

- Acquired throughout Summer of 1973 

- Six 5x20-mile segments in Indiana and Illinois 

- Cloud cover substantially reduced the number of useable acquisitions, 

Fayette, III., Segment had the most, useable acquisitions 
(Early and Late June, Mid July, and Late August) 

- Color composite images recently produced on the LACIE Production Film 
Converter (PFC) 

- LACIE-sized (5x6-mi) segments recently extracted from the 5x20-mile segments 

Field Identification Information (^^25% of Area) 

- 20 Quarter sections by periodic ground observations 

» 

- 20 Full sections by interpretation of multidate aerial photographs 
Agronomic Data 

- Periodic (18-day) ground observations of plant height, ground cover, 
crop maturity stage, and condition 

- Field size, row width, and row direction 



TABLE 7.4 


LANDSAT DATA ANALYZED 


• Pixels Within Identified Fields 

“• Boundary or edge pixels were excluded to eliminate this source 
OF confusion from initial studies 

- Crop label and field number were appended to each pixel in 
Fayette segment 

- Multi temporal data set assembled for Fayette segment 

• Landsat Mean Vector From Each of These Field Centers 


~ Mean vectors used in CITARS were analyzed 
- All segments and all useable dates 



TABLE 7.5 


LACIE PHASE 3 BLIND SITE DATA 


Eleven Segments Prepared for Analysis 

(Most included in JSC's list of candidate segments for multicrop analyses,) 

- Screened (visually and by computer algorithm) 

- XSTAR preprocessed (sun angle and haze corrections) 

- Extracted 897 Points (5x5 gridj includes 209 dots as subset) 

- Merged with ground truth data from JSC tapes 

(Labeled as ei'ther a single specific cover code or a mixture pixel.) 


Segments^ With Corn and Soybeans Percentages (From 897 Points) 


State 

Segment 

% Corn 

% Soybeans 

State 

Segment 

% Corn % 

ND 

1619 

2.3 (21)* 

0 

SD 

1669 

1.2 (11) 


16A0 

3.2 (29) 

0 

• 

1686. 

10.9 (98) 


1663 

O.A (A) 

5. A (A8) 


1699 

2.6 (23) 


1903. 

2.6 (23) 

0 


1802 

15.5 (139) 


1927 

5.2 (A7) 

0 


1805. 

7.6 (68) 

MINN 

1515. 

A. 3 (39) 

3.6 (32) 





Number of points. 


Soybeans 

0 


0 


CD O O 





TABLE 7.6 

PLANS 


• Apply XSTAR Preprocessing and Tasselled-Cap Transformation to Field-Center Data 
From All Segments and Analyze 

(Calibrations of the early Landsat-1 CITARS data and the later 
Landsat-1 and Landsat-2 data have been compared under the other 
ERIM Multicrop Task.) 

• Perform Quantitative Analysis of Periodic Observations From All Segments 

- Relative crop calendars 

- Relative crop growth and condition 

- Relate to spectral signatures 

• Initiate Analysis of Separability in Larger Contiguous Ground Areas^ Including 
Boundary Pixels 
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CITARS DATA SITES 


O 




ERTS-1 ' C'snS / 
passes: 0/P 


X) I K\\N\1 / 

H/H / 



study Area Counties: 
Indi ana 


N ' \ 'One segment : 

'/a X 32 km 
25,856 hectares 
(64,640 acres) 


One sectlon:- 
256 hectares 
(640 acres) 


n 






m 



ERTS-l 
overl ap 


1 1 H nol s 


Data Acquisition Periods: 

0-5/21-25/73 IV - 8/01-05/73 



1. 

Huntington 

4. 

Ll v1 ngs ton 


2. 

Shelby 

5. 

F aye tte 


3. 

Whi te 

6. 

Lee 


* 

Ground 

Truth : 




ASCS 

— 20 quarter secti 


I - 6/08-12/73 
n - 6/26-30/73 
in - 7/14-18/73 


V - 8/19-23/73 
VI - 9/06-10/73 
VII - 9/24-28/73 


Photointerpretation — 20 sections (black) each ERTS-1 pass 




BAND 6 Ce^lO) 


FIGURE 7.2 


MULTITEFIPORAL CLUSTER ANALYSIS - FAYETTE COUNTY 

Pixels Labeled Corn; Four Dates, 16 Channels 

Pattern of Corn Clusters on First Two Dates 
(338 Pixels) 




BAND 5 (6/10) 


BAND 5 (6/29) 


10 June 1973 
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FIGURE 7.3 


MULTITEMPORAL CLUSTER ANALYSIS - FAYETTE COUNTY 

Pixels Labeled Corn; Four Dates^ 16 Channels 

Pattern of Corn Clusters on Last Two Dates 
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FIGURE 7,4 


MULTITEMPORAL CLUSTER ANALYSIS " FAYETTE COUNTY 
Pixels Labeled Soybeans on 21 August; Four Dates^ 16 Channels 


Pattern of Soybeans Clusters on First Two Dates 

(508 Pixels) 
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FIGURE 7.5 


HULTITEUPORAL CLUSTER ANALYSIS - FAYETTE COUNTY 
Pixels Labeled Soybeans on 21 Augustj Four Dates . 16 Channels 

Pattern of Soybeans Clusters on Last Two Dates 
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FIGURE 7.6 


UNITEffORAL CLUSTER ANALYSES - FAYETTE COUNTY 

Results Using Pixels Labeled Corn or Soybeans 
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UNITEMPORAL CLUSTER ANALYSES - FAYETTE C0U|TY 

Results Using Pixels Labeled Corn or SoYBEAPis 
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FIGURE 7.8 


UNITEMPORAL CLUSTER ANALYSIS - FAYETTE COUNTY 
Result$ Including Pixels From Other Crops, 17 July 1973 Acquisition 
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UNITEMPORAL CLUSTER ANALYSIS - FAYETTE COUNTY 

Results Using All Labeled Pixels/ 21 August 1973 Acquisition 



tojER., Name 
C.QRR Soybeans 


Number of Pixels; 
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Composition (%): 



Label Crop 

66 
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Other Crops 

35 

13 

% of Major Crops; 

72 
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(Number of pixels) 
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FIGURE 


MULTITEFIPORAL CLUSTER A 

Results Including Pixels from 0th 

Patterns of Clusters 
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FIGURE 7.11 


PIULTITEMPORAL CLUSTER ANALYSIS - FAYETTE COUNTY 

Three Dates? 29 June^ 17 July^ 21 August 1973 

Results Including Pixels From Other Crops (continued) 
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CROP CODES: 
CODE CROP 


C 

S 

T 

I 

H 

W 

X 


CORN 

SOYBEANS 

TREES 

CITY 

HAY 

WHEAT 
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10 


5 (Aug) 


80 


Pattern of Clusters on 21 August 1973 


Cluster Statistics 


Number of Pixels? 

Composition (%)? 
Label Crop 
Other Crops 
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o p 

St5 
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91 

92 

91 

92 

9 

8 

'84 
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« 1380 pixels total- clusters shown contain 92%of tota 
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TASK 8 

MULTICROP LABELING AIDS 
(R. C. Cicone, P. F. Lambeck, and R. Balon) 

8.1 INTRODUCTION 

The accuracy of operational large area crop inventory systems 
modeled after LACIE depends critically on the correctness of crop 
labels generated by Analyst Interpreters. Task 8 is being conducted 
as a partial response to a request from the multicrop inventory plan- 
ning committee for support from the SR&T community in adapting LACIE 
technology to a multicrop environment. The critical issue addressed 
is that of analyst labeling. 

8.2 OBJECTIVES 

Our purpose through Task 8 is to analyze the methods of presenta- 
tion of Landsat data for purposes of human interpretation to assess how 
well they convey the information relevant to crop discrimination. We 
endeavor to develop new data presentation techniques in the form of 
false color image products as well as graphic displays of spectral infor- 
mation, which stand to expedite correct labeling of crops. Task 8 is 
divided into two subtasks. 

Subtask One has the objective of evaluating and making recommendations 
for false color image products suited for multicrop labeling purposes. 

Subtask Ttto has the objective of developing and evaluating spectral 
aids to assist analyst interpretation of multicrop classes . 

8.3 APPROACH 

Points of approach for this task are layed out in Table 8.1. 

8.4 SUMMARY OF PROGRESS FOR THIS QUARTER 

The major accomplishments of this quarter are listed in Table 8,2. 

Four issues were addressed and will be discussed in the following 
section. 
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8.5 DESCRIPTION OF PROGRESS 

The first issue addressed this quarter was an analysis of the 
sensitivity of image products to the information content of Landsat 
data. A method for graphic display of image sensitivity by analytic 
means was developed during the quarter. 

The technical issue under scrutiny is whether the information con- 
tent in a false color image as viewed by an AI is consistent with the 
information content of that same scene in digital form. Table 8.3 
describes the technique for display of image sensitivity. To begin 
with, a distance in color space of five counts was empirically deter- 
mined (using false color image products) to represent a perceptible 
difference in color of objects. A grid of sample points was then 
placed in the brightness/greenness plane of the Landsat digital data. 
This grid was placed at a six count interval and spanned the range of 
the color production capabilities of the production film converter 
(i.e., any point outside this range Would be saturated in color). One 
can picture circles centered at each grid point with a radius of three 
counts and tangent to adjacent circles. These three count circles ^?ere 
chosen so as to conservatively contain the noise characteristics of the 
Landsat data about the grid points. The question addressed, then, was 
CO determine whether sufficient color resolution was provided, in light 
of a five count color resolution capability, to preserve information 
in the Landsat digital data at the level indicated by the circles. 

Figure 8.1 illustrates a typical result using Product 1 color 
mapping procedures. Each ellipse represents visual sensitivity to 
color changes about the center. Resolution of the data into color is 
poor along the direction of elongation. Smaller ellipses mean the data 
space is better represented by separable colors. Large ellipses indi- 
cate that color varies slowly. It is certainly clear that three count 
tangent circles are not preserved in this mapping. Figure 8.1 Indicates 
poor resolution of the data into color and 'distortion of the information 
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in certain portions of the space. Figure 8.2 overlays Figure 8.1 to 
show where principal crop distributions corresponding to this acqui- 
sition fall. 

We have found that the characteristic pattern of sensitivity results 
from the exposure control linkage of the Production Film Converter, 
which relates input level exponentially to film transmission. This 
results in decreasing sensitivity toward less than the highest input 
levels, or, equivalently, the highest Landsat channel values. Overall 
image sensitivity is very much a function of the bias and scale factors 
selected for data transformation before image generation. It was dis- 
covered that the image sensitivity resulting from use of transformation 
factors appropriate to an acquisition displaying the full natural range 
of agricultural data, was inadequate for resolving differences of less 
than 7 data counts. 

A hypothetical PFC calibrated to relate input level to film trans- 
mission linearly (rather than exponentially) was investigated. Image 
sensitivity was found to be lower. Figure 8.3 shows the sensitivity 
display produced in this experiment. 

A second study conducted this quarter was a probability of mis- 
classification analysis of crop distributions processed through various 
image generation procedures. This was done with the idea of determining 
if and to what extent loss of information occurred in the data-to-false- 
color-image transformation. The crop distributions used were wheat and 
non-wheat clusters for the Kansas sample segment 1154, acquisition of 
day 124, 1976 (Figure 8.2). This acquisition was chosen for the good 
separation between these txTO classes. 

We found that dropping a channel of spectral information (Landsat 
Channel 3) caused a sizeable increase in expected misclassif ication 
betxjeen the crop distributions. Table 8.4 presents details. This has 
been recognized empirically and the AI is routinely provided a product 
which includes Channel 3 (Product 2), Hox-rever, the price is that this 
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image does not contain Channel 1 information. We expect that much of the 
information of Landsat data lies in the interrelationship of the four 
channels- considered simultaneously. For this reason, interpreting two 
images with subsets of bands is not the same as interpreting one image 
with all the information. Data feature extraction transformations exist 
which reduce the dimensionality of Landsat data without significant loss 
of information. The linear Kauth transformation is a prime example. 

There is impetus for investigating ways of converting Landsat features 
to a single color image with full retention of information. 

The two crop distributions were converted into corresponding 
color distributions, simulating the Product One algorithm and proposed 
variations as well. Tables 8.5 and 8.6 present details. In each case 
an increase of confusion between the distributions was observed x^hich 
was of comparable seriousness to that due to dropping a band. This 
confusion is attributable to incomplete resolution of the information 
in the data into color. We have examined PFC products for which color 
coordinates of field means have been computed. From this we judge the 
human eye cannot xjell distinguish colors separated by 5 color space 
units or less (x^e use the 1976 CIE Uniform Color Space), For the present 
VTB use 5 units as our standard of color resolution. If an image product 
takes meaningful data variation into perceptually indistinguishable 
colors then information is lost. Our sensitivity analysis presented 
above suggested this was happening for a typical scene, and our classi- 
fication accuracy study tends to confirm that conclusion. 

Our third activity during the quarter related to facilitating the 
transition to multicrop interpretation and labeling. The variation of 
corn and soybeans over time xxras plotted in color space coordinates. 

The data used \<ias from CITARS acquisitions of Fayette County. Figure 8.4 
shoxre color trajectories of the crops in the principal components of 
color variation. 
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A fourth activity, in response to reports of apparent problems in 
applying the SCREEN and XSTAR algorithms to CITARS data, determined that 
the apparent problems observed by ERIM and by UCB were caused in each 
case by errors in processing. Our analysis indicates that the Landsat-1 
to Landsat-2 data conversion previously specified by ERIM is actually 
very close to the optimum conversion for CITARS data. SCREEN and XSTAR 
performed normally on this data. 

8.6 CONCLUSIONS AND RECOMIENDATIONS 

Specific conclusions are listed in Table 8.7. We recommend investi- 
gation of new data-to-image processing to minimize loss of information 
content . 

8.7 PLANS 

Plans for the coming quarter are detailed in Table 8.8. 
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Table 8.1 

APPROACH 


SUBTASK 1 

® Use Principles of Color Science to Model the Color Production Characteristics 
OF the PFC 

® Analyze Techniques for Producing False Color Imagery from Landsat Data in 
Terms of: 

— Spectral information content 

00 

— Variation of color sensitivity in image data space 
* 

©Evaluate Trends of Color Variation for Corn and Soybeans as Function of Date 
©Examine Alternative Single and Multi-Date Image Products 

SUBTASK 2 

® Evaluate Existing Multicrop Data Sets to Examine; 

— Effects of haze and sun angle compensation techniques 
— Develop spectral aids parallel to those currently employed in Lacie 

©Examine Utility of Haze and Sun Angle-corrected Image Products 


iit atud 



Table 8.2 


SUmARY OF RESULTS 

® A Method for Graphic Display of Sensitivity of Film Products Throughout 
Data Space was Developed and Used to Evaluate Image Sensitivity, 

® An Analysis of Image Products was Performed in Terms of Probability of 
Misclassification of Crop Distributions in a Color Space that is 
Perceptually Uniform with Respect to the Human Eye, 

® Preparation of CITARS Data Set, 

o Growing Season Color Trends for Corn and Soybean Were Examined Using 
CITARS *Data. 

o Performance of the SCREEN and XSTAR Algorithms on CITARS Data Was Verified. 
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Table 8.3 

IMAGE SENSITIVITY DISPLAY PROCEDURE 


• Sample Points in Data Plane at a Spacing of 6 Counts 
•Transform Them to Uniform Color Space 

• Construct a Symmetric Normal Distribution About Each Point With a 
Standard Deviation of 5 Units 

• Analytically Transform Each Distribution Back to the Data Plane 
AND Display 
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Table 8.4 


INFORMATION LOSS DUE TO USING A SUBSET OF FEATURES 
DURING BIOWINDOW 2 (S/S 115A) 


Percent Misclassification 

Variables Wheat Other 


Bands h 2, 3 and 4 

2.. 8 

5.4 

Bands Ij 1 , and 4 

12.2 

4.8 

Tasselled Cap Brightness/Green 

10.4 

4.6 
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Table 8 . 5 


INFORMATION LOSS DUE TO COLOR RESOLUTION OF THE HUMAN EYE 




Percent Misclassification 



Wheat 

Other 

Product 1 

a = 0 

15.2 

5.0 

Product 1 

a - 5 

18.2 

8. A 
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Table 8.6 


INFORMATION LOSS DUE TO FILM PRODUCT 
Biowindow 2 (S/S 1154) 

Percent Misclassification 

Product and Processing Wheat Other 


• No Preprocessing 


Product 1 

18.2 

8.4 

Kaneko Product 

18.4 

17.4 

Kraus Product 

26.4 

32.8 

Double Band 4 

Kaneko Product 

17.2 

6.2 

Hocutt Product 

17.2 

6.6 

Kaneko' s Processing 

Kaneko Product 

20.6 

8.6 

Hocutt Product 

18.6 

6.4 

Kraus Product 

21.2 

9.6 


m 
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Table 8.7 


CONCLUSIONS 


• With Band 4 Doubling, Kaneko, Hocutt and Product I's Retain Comparable 
Information and Display Similar Color Sensitivity 

• Kaneko Preprocessing Resulted in Slight Deterioration of Expected Classi- 
fication Using Alternative Productsj Kraus Product Lags in the Comparison 

• No, P reprocessing Resulted in the Poorest Expected Classification 

f 

• With Either Exponential or Linear Film Transmission Control, Color Resolution 
of the Eye Masks Available Resolution in the Data in the Region of Data 
Concentration 

•Bias and Scaling Techniques Utilizing Bands 1, 2, and A As Inputs to the PFC 
Color Guns Reduces Information Content of the Data and Does Not Preserve 
Distance' Relationships in the Data 
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Table 8.8 


PLANS 


• Prepare an Interim Technical Report 

• Examine Image-to- Image Color Characteristics of Crops 

• Develop and Evaluate Techniques for Data-to-Image Happing 
Which Will Retain Information Content of Imagery 

• Consider Hulti temporal Image Products 

• Begin Development of AI Spectral Aids Based on Current Results 
From Crop Discriminability Task 

• Investigate Application of XSTAR Haze Correction to Address 
Localized Atmospheric Variations Within the Boundaries of 
LACIE Segments 



Figure 8.1 







Figure 8,2 


Data Location for S/S 1154 acquisition 
in Bio 2 



OTHER 
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Figure 8.3 

SEGMENT 1154 BIO 2 
PRODUCT ONE WITH LINEAR PFC CONTROL 
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